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I. INTRODUCTION 

As in any conductor, the potential due to a source dis- 
tribution is screened within an electrolyte solution. The 
screening charge carriers, or counterions, in these classi- 
cal systems are ionized atoms or molecules. The quali- 
tative result of the screening is typically only to reduce 
the distance out to which the potential from a source 
distribution is significant. However, experiments have 
shown that multivalent counterions can allow for more 
exotic behavior For example, over-efficient screening 
due to multivalent counterions can sometimes result in 
the charge reversal of a molecule or substrate. In a re- 
lated effect, multivalent counterions have been observed 
to mediate attractions between similarly charged macro- 
molecules. Study of these counterintuitive effects is perti- 
nent to the understanding of many biological and colloid 
systems, in which molecular conformations and interac- 
tions are often sensitively controlled by counterion con- 
centration. 

The Poisson-Boltzmann (P-B) mean-field theory pro- 
vides one of the simplest descriptions of electrolyte re- 
sponse and often allows for accurate modelling of these 
systems. However, it has been proven for mirror- 
symmetric geometries that like-charge attraction is im- 
possible for a pair of identical molecules within this the- 
ory This has bolstered the general opinion that the 
physical mechanisms behind like-charge attraction and 
overcharging cannot be understood within a mean-field 
description of electrolytes. Most recent works have thus 
attempted to understand these observations through a 
consideration of fluctuation corrections to the mean-field 
response Certainly, for homogeneous systems, 

it seems clear that fluctuations provide the only pos- 
sible mechanism for bringing about the charge correla- 
tions necessary for these effects. This conclusion does 
not immediately follow for heterogeneous systems, how- 
ever, since static modulation allows for the introduction 
of some inherent charge correlation. Indeed, studies have 
shown that heterogeneity can play an important role in 
these systems. For example, numerical simulations have 
shown that modulation of molecular geometry can sig- 



nificantly affect interactions [15|. In addition, it has 
been shown perturbatively that within the P-B theory, 
charge condensation onto a surface in an electrolyte so- 
lution may be increased when the surftice charge contains 
certain limiting forms of modulation [16]. This extra con- 
densation can in turn lead to a decrease in the osmotic 
pressure between two surfaces, demonstrating that mean- 
field effects can at least help to facilitate both overcharg- 
ing and like-charge attraction [l7j . 

It is the purpose of the present paper to further explore 
the effects of inherent, quenched modulations within elec- 
trolyte solutions. Fluctuation effects are neglected and 
the electrolyte response is considered within the P-B for- 
malism. Significantly, it is shown that inherent modu- 
lations alone, can, in fact, allow for both like-charge at- 
traction and overcharging mechanisms. The key to like- 
charge attraction within the P-B theory is the break- 
ing of mirror symmetry, which was assumed to hold in 
Ref. [H. A simple example is given which demonstrates 
this point. Heterogeneity- induced condensation is stud- 
ied using a perturbation approach similar in spirit to that 
considered in Ref. fl^ . The present approach consists of 
expanding about solutions to the linearized P-B equa- 
tion and allows for the consideration of additional charge 
modulation limits, boundary surface modulations, and 
the effects of having more than one species of screening 
ion. Various examples are considered in order to provide 
a brief survey of these different set-ups. Of particular 
interest is the counterion concentration dependence ex- 
hibited by systems containing multiple species of ion. In 
principle, this dependence should often allow one to de- 
termine whether modulation is a significant contributor 
to overcharging in any particular experimental observa- 
tion. 

An outline of this paper is as follows. In section II the 
P-B theory is reviewed and solutions to the linearized P- 
B, or Debye-Hiickel 18], equation are presented. Charge 
and surface modulations are considered in sections HI 
and IV, respectively, while section V contains concluding 
remarks. 



II. THE POISSON-BOLTZMANN EQUATION 
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If one assumes that the mobile charge distribution in 
an electrolyte solution is given by a sum over ion species 
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of Boltzmann factors multiplying bulk charge densities, 
then the Poisson equation reads 



i 

« 72^-47rp,. (1) 

Here, is the electrostatic potential, the qi and rii are the 
charge and bulk densities of screening species i, kT is the 
thermal energy, pg is the source distribution, and 7 = 
is the inverse Debye length [l^ . The first equation above 
is the full P-B equation. The linearized approximation is 
written out in the second line. Solutions to these equa- 
tions provide mean-field approximations to the potentials 
of physical systems. 

It is possible to derive the P-B equation in a rigorous 
manner which clarifies the physical approximations and 
assumptions required to arrive at this expression. In the 
weak potential limit, when the linearized form applies, it 
has been shown that the equation is exact, while errors 
result at each non- linear order [2^[2l|. However, recent 
simulations have shown that a partial cancellation of two 
neglected effects, counterion correlations and finite ion 
sizes, can result in approximate agreement between the 
P-B solutions and experiments |22| . Scattering experi- 
ments have also shown that the full non-linear solutions 
are in reasonable agreement with observed charge distri- 
butions when only a single screening ion species is present 
[23| , though some extra divalent condensation was ob- 
served when both monovalent and divalent ion species 
were present [l^]. These results justify use of the non- 
linear equation, but the solutions should be considered 
qualitatively accurate only outside the innermost. Stern 
layer surrounding highly charged distributions. 

Solutions to the inhomogeneous, linearized P-B equa- 
tion will be needed throughout this paper. The Poisson 
sum rule provides a convenient method to determine the 
potential of a lattice of point charges, from which one 
can obtain the potentials of more general distributions 
through superposition. For an arbitrary individual par- 
ticle potential function /, the rule states 

^ /(r - R) = i / ^ /(r - R) cxp(zG • R)d"i?. (2) 

R G 

Here, n is the dimension of the lattice, V is the n- 
dimensional volume per unit cell of the direct lattice, the 
vectors R are the direct lattice vectors of the distribu- 
tion, and the vectors G are the reciprocal lattice vectors 
of the distribution 25] . Shifting the origin of integration 
in Eq. ^ to the projected position of the observation 
point r immediately gives formal integral representations 
for the Fourier coefficients of the potential. 

To obtain the potential from a lattice of Coulomb point 
charges which are linearly screened, one may plug in the 
Yukawa potential function / — (7exp[— 7r]/r, which is the 
solution to the linearized P-B equation with point charge 



source q. The resulting expressions for the potentials of 
one, two, and three-dimensional Yukawa lattices are [2^ 

<I>i,y = 2A^i^o((|G|2-K7')i/2d)exp[zG-r], (3) 
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where A, cr, and p are the average charge densities of 
the lattices. For large arguments, the modified Bessel 
function Kq above may be expanded as [l^l 



V 2z ^-^ 



(6) 



It follows that each component will be exponentially 
damped with both the frequency and the distance from 
the distribution for one and two-dimensional systems. To 
relate these expressions to their unscreened analogs, one 
need only replace (jGp -|- 7^)1/2 by |G| in the non-zero 
frequency components and then also replace the zero fre- 
quency component by the appropriate continuous charge 
distribution expression. This is given by — 2Alogr for a 
linear distribution, for example. 



III. SOURCE CHARGE MODULATION 

Consider now the charge modulated system depicted 
in Fig. [T] Two model polyelectrolytes are shown, each 
of which consists of alternating point particles of charge 
-f 2 and —1 fixed on a rigid backbone. The period of the 
modulation is taken to be 1. Note that each polyelec- 
trolyte has a net charge of -1-1 per period. If the two 
polyelectrolytes are shifted with respect to one another 
by half of a period, as shown in the figure, then they 
will attract one another when their separation is small, 
demonstrating that mean-field like-charge attraction is 
possible. 

One can approximate the interaction energy of this sys- 
tem in the Debye-Hiickel limit by superposing the poten- 
tial for an infinite line of —1 charges with that for a line 
of +2 charges, each obtained by plugging into Eq. ([3]), 
and truncating the series after two terms. The resulting 
approximate interaction energy per period of the second 
polyelectrolyte is given by 



E{d) « 2K0M - 36Ko[v/(27r)2+72d], 



(7) 



where d is the separation distance between the two poly- 
electrolytes. A plot of this function is shown in Fig. [2] 
for 7 — 3.0. Attraction is indeed observed at small 
separations while repulsion is observed at large separa- 
tions. This form of attraction is not salt mediated, but 
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FIG. 1: Portions of two long, identical polyelectrolytes are 
shown. At large distances they repel but at short distances 
they attract. 
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FIG. 2; A plot of the approximate interaction energy versus 
separation for the model polyelectrolytes. 



as 7 is increased, the repulsive effect at large distances 
is reduced. This trivial example demonstrates that like- 
charged molecules can be attracted to one another even 
when the electrolyte is considered at the mean-field level. 
Note, however, that if the model polyelectrolytes are ar- 
ranged in a mirror-symmetric orientation, as was explic- 
itly assumed in Ref. [2] , then the coefficient of the second 
term in Eq. ([7]) would also be positive. This is consistent 
with the conclusion that identical molecules cannot be at- 
tracted to one another, within the P-B formalism, when 
they are arranged in a mirror-symmetric orientation. 

We now turn to a consideration of charge modulation- 
induced condensation. As mentioned above, this was 
previously studied in Refs. [l^ and These previ- 

ous works considered solutions to the full non-linear P-B 
equation for systems containing a planar source distri- 
bution and only a single species of screening ion. Solu- 
tions were obtained up to second order in the parameters 
e(q)/CT, for all g ^ 0. Here, a was the average charge 
density of the planar source distribution and the e{q} 
were the amplitudes of the oscillatory components of the 
source distribution's Fourier decomposition. Thus, these 
solutions were valid in the limit of relatively weak modu- 
lation relative to net surface charge. Interestingly, these 



solutions indicated that each non-zero Fourier component 
independently led to an increase of screening charge con- 
densation at second order in perturbation theory. This 
universal result was also shown to hold for long wave- 
length modulations. Here it is shown that a similar con- 
clusion also holds true in the general weak potential limit. 

For simplicity, we consider first a system which con- 
tains only one species of counterion. The counterions are 
taken to have charge q and bulk density ng, so that the 
P-B equation reads 



V = —AnquQe 



-q4>/kT 



(8) 
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Here, ps is a perturbing planar charge distribution. The 
first term in the above Boltzmann factor expansion cor- 
responds to the uniform background charge density. This 
does not effect an electric field and so may be dropped 
for now. We assume that the potential can be expanded 
in a series a.s <p = (pi + (j)2 + where the term (/j^ is 
of order p^. Collecting like-powers of Ps then gives the 
system of equations 
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The solution to the first equation above may be obtained 
from our earlier expressions if the distribution ps is ex- 
panded in a Fourier series as 
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Ac exp[iG • r]. 



(10) 



Note that while extra, oscillatory condensation results 
at this first order, it averages to zero when integrated 
over the modulation directions. The source term for (^2 
is proportional to the square of the first order potential, 
however, and a non-zero condensation will result at this 
order. To obtain the net condensation, one need only 
consider the average charge density above the surface. 
The equation for <j)2 is therefore averaged over x and y 
to obtain 
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an ordinary differential equation which may be solved 
trivially. The solution is 



(12) 



X {exp[~7z] 
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exp[-2\/G2-K 72 2]}. 
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A plot demonstrating the resulting z dependence of a 
typical term of the second order average charge density is 
shown in Fig. |3l The total averaged second order counte- 
rion charge density at the surface is given by the following 
sum over modulation components: 



P^(2 = 0) = - 
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(13) 
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This demonstrates that indeed there is extra condensa- 
tion proportional to the square of the amplitude for each 
mode of the modulation, consistent with the long wave- 
length and small modulation results of Refs. 



161 and [U 



The present analysis extends the conclusion that univer- 
sal modulation-induced condensation is expected at sec- 
ond order to all relative values and wavelengths of modu- 
lation. In particular, the result applies to systems which 
have little or no net charge but large modulations. In 
systems such as these, the modulation-induced conden- 
sation could easily result in overcharging. This demon- 
strates that mean-field overcharging is possible in certain 
limits. 

A single component electrolyte is somewhat non- 
physical in that the solution contains charges of only one 
sign and thus cannot be net neutral. Accordingly, it is 
of interest to consider systems which contain multiple 
species of ions. This would be difficult when considering 
solutions to the full P-B equation but the present solu- 
tion technique is readily applied to such systems. In a 
general, neutral electrolyte solution, the potential satis- 
fies the approximate equation 



(14) 



Here, the coefficients 7^ and C are now given by sums 
over ions species as 
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(15) 
(16) 



screening charges of both signs and equal magnitudes. 
For this system, the quadratic terms in the expansion of 
the Boltzmann factors will cancel and no second order 
potential or condensation will result. This implies that 
any extra condensation would have to occur at higher 
order for such systems and the resulting effect should be 
much weaker. The quadratic term in the expansion is re- 
tained, however, for a net neutral system which contains 
charges of different sign and different valence. For exam- 
ple, if the positive counterions have charge qi with bulk 
density hq and the negative counterions have charge ~q2 
with bulk density ^?^o, then the resulting second order 
coefficient is given by 



_ 27rgino , 2 



(18) 



The result is that extra charge condensation still occurs 
at the lowest non-linear order for a net neutral system, 
provided the free charges have different valences. Note 
that the net sign of the resulting modulation-induced 
condensation does not depend on Aq. Extra condensa- 
tion occurs for both species of ion, and the sign of the 
net condensation depends only on which species domi- 
nates. This, in turn, is determined by the sign of the 
quadratic term in the expansion of the P-B equation; if 
C is large and positive, the resulting condensation will be 
large and positive, and vice versa. This obviously holds 
for any distribution of counterion species. Further com- 
ments on this point appear in the discussion section of 
the paper. 

One and three-dimensional charge modulated distribu- 
tions may be analyzed in a similar manner. For a three- 
dimensional lattice, the concept of condensation is not 
well-defined. In the one-dimensional case, however, an 
analogous condensation to that discussed above for sur- 
face distributions also occurs. Each modulation compo- 
nent leads to additional condensation at second order, the 
sign of which depends only on the sign of the quadratic 
term in the expansion of the P-B equation. A considera- 
tion of similar surface modulation effects is contained in 
the following section. 



Charge neutrality requires "^qiUi = 0, if the charge of 
the polyelectrolyte itself is neglected, for simplicity. Solv- 
ing the resulting system of equations as in the above gives 
the following for the average electrolyte charge density at 
the surface: up to second order. 
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where po is the charge density which would appear for a 
uniform distribution of charge. 

It is interesting to consider some particular examples. 
Consider first a net neutral electrolyte system containing 



IV. SURFACE MODULATION 

In the above, we saw how charge modulation can al- 
low for oscillatory potential components that lead to a 
net charge condensation at non-linear orders. A modu- 
lated surface can also allow for the introduction of po- 
tential modulation and similar condensation effects must 
naturally result. It is worthwhile to consider some exam- 
ples, and two are presented here which are susceptible to 
perturbation expansions. These examples are of interest, 
not because they accurately model physical systems, but 
because they demonstrate the mechanisms involved for 
different boundary conditions and geometries. 
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FIG. 3: A plot of the one component fluid's second order 
charge condensation versus the distance from the plane. 
Here G = 2'^ = 2/Id ■ The condensation is particularly local- 
ized and then screened at larger distances. The localization 
increases exponentially with the magnitude of G. The total 
integrated charge from z = to oo is zero for each component, 
as expected. 



A constant potential boundary condition problem will 
be considered first, where the surface takes the form 



iG-i 



(19) 



Here /? is a small constant and the Ac are assumed to 
be 0(/3°). At hnear order, the system of equations to be 
solved is 

'(V^^t')*^! =0 \iz>Z{x,y), 
4>i = l iiz = Z{x,y), (20) 

— > as z oo, 

and the term (jji is assumed expandable as 

oo 

'^1 =$]<^i,fe/3'. (21) 

fc=0 

Each of the (pi^k satisfy the Debye-Hiickel equation indi- 
vidually and may be determined iteratively by expanding 
the boundary condition as 



1 = 0i(z = Z) 

EfSr^ A „»GT\n^" ^^^^ 



n=0 G 



i! dz" 



(22) 



Equating coefficients at order on both sides of Eq. (|22|) 
generates the boundary conditions for the (pi^k at the 
simple surface z = 0. Solving up to second order in /3 
gives 
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The last term has once again been averaged over the mod- 
ulation directions to simplify the calculation. As in the 
charge modulation discussion, each modulation compo- 
nent is observed to independently result in extra second 
order condensation. Here, however, the surface modu- 
lations generate the extra charge condensation at linear 
order in the P-B theory. Since the effect takes place at 
linear order, the sign of the condensation does not de- 
pend on the sign and distribution of the counterions, but 
instead is always opposite to that of the surface potential. 

The first order solution above can easily be plugged 
back into the non-linear P-B equation to obtain the next 
order term in the solution. Although this won't be done 
here, we note that the non-linear terms in the P-B ex- 
pansion can result in extra condensation at the same or- 
der in the potential and in /3 as that obtained above. 
These terms, however, are sensitive to the sign of In 
particular, if C = 0, the linearized P-B solution will be 
correct to lowest order. Alternatively, the higher order 
solutions can result in charging of the wrong sign, or neg- 
ative screening, if ^ is chosen appropriately. 

We now consider a nearly cylindrical geometry with 
boundary conditions which are more applicable to bio- 
logical molecules. The set-up consists of a line of charge 
of unit linear charge density surrounded by a dielectric 
material of shape 



i?= 1 + /3cos(Gz). 



(26) 



Outside the dielectric near-cylinder the potential is as- 
sumed to satisfy the linearized P-B equation. The fol- 
lowing system of equations specifies the unique solution 
for the potential: 



f(V2 



-f^)(t) = if r>i?(z), 

-^p if r<i?(z), 

0-^-0 as r — > oo, 

4>i ^ 4>o at r = i?(z), 

T>i ■ n = Do • fi dX r — R{z). 



(27) 



Here, the subscripts i and o stand for inside and outside, 
respectively, e is the dielectric coefficient, D is the elec- 
tric displacement, n is the outward normal to the near 
cylinder's surface, and p = <5(r)/27rr. 

The boundary conditions may be expressed at the sur- 
face r = 1 in a similar manner to that presented in 
Eq. ((22|) . For example, the boundary condition relat- 
ing the equality of the normal component of the electric 
displacement inside and outside the dielectric-solution in- 
terface takes the form 
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where the normal is given exphcitly by 

r l3Gsm{Gz)z 



A eff 



n 



(29) 



Up to second order, the boundary conditions generate 
terms of the form 



2A 

A logr + l3BIo{Gr) cos(Gz) 



+I3^\C + DIo{2Gr) cos(2Gz) 



(30) 



and 



(bo = EKohr) + pFKoiiG^ + jY^^r) cos(Gz) 

+/32|Gifo(7r) + ^f^o((4G2 +72)1/2^) cos(2Gz)|, 

(31) 

where Iq is the modified Bessel function of the first kind. 
The boundary condition equations further enable one to 
obtain a system of linear equations which specify the 
unique solution for the coefficients A through H. The 
algebra is most easily carried out in a software package 
such as Mathematica. The resulting expressions are easy 
to obtain but are too lengthy to report here. Once the 
coefficients for the potential have been determined, one 
can read off the effective charge of the near-cylinder, de- 
fined here to be {E + /S'^G) /2. This is proportional to the 
coefficient of Ko{'^r) above, a natural choice given that 
the potential due a screened line of charge is given by 
2XKoi-/r), from Eq. ©. 

We now examine how the resulting effective charge 
varies with the wavelength of the surface modulation. 
At long wavelength the surface curvature is low and the 
cylinder has a slowly varying local radius. The effective 
charge at long wavelength is thus expected to be that 
which one would obtain through an appropriate aver- 
age over radii of the effective charge for an unmodulated 
cylinder. The appropriate average is 



du 



760(1 -f /3 cos 27ru)iiri [7(1 -f /? cos 27ru)] 



(32) 



This integral is over a convex function, which, by Jenson's 
inequality, implies that the modulation will result in an 
increase rather than a decrease in the effective charge 
of the cylinder. A plot of the effective charge, valid to 
second order in /?, versus wavelength is shown in Fig. 2] 
The effective charge does approach the radii-averaged un- 
modulated value at long wavelengths, as expected. Note 
also that the effective charge decreases as the modula- 
tion wavelength decreases, and below wavelength values 
oi Iq ~ 2, the resulting condensation works to reduce 
the net effective charge of the cylinder. However, the 
perturbation expansion converges most quickly when the 
modulation wavelength is large. Indeed, the expansion 
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FIG. 4: A plot of the effective charge density versus modu- 
lation wavelength for the near-cylinder. The thin line is the 
unmodulated charge density, the dashed line is the modulated 
effective charge valid to second order in /3, and the thick line is 
the radii-averaged value discussed in the text. The parameter 
values were set to /? = 0.1, ei = 1, eo = 5, and 7 = 5. 



will break down when /3G > 1 owing to the fact that the 
/? power series of n in Eq. (j29p will diverge in this limit. 
Thus, convergence concerns indicate that the plot should 
not be considered quantitatively accurate at small values 
oi Iq. a decrease in the effective charge is also observed 
at short modulation wavelengths for very small /3, how- 
ever, where the second order solution is likely to be very 
accurate. It is thus reasonable to expect the observed 
decrease in effective charge at small modulation wave- 
lengths to be qualitatively accurate at larger f3 values as 
well. 

These two examples demonstrate that the sign of sur- 
face modulation-induced condensation can be sensitive 
to both counterion valence and modulation wavelength. 
A modulation wavelength dependence was observed only 
in the second example but it also holds true for some 
other systems, linear geometries with constant potential 
boundary conditions providing one example. We have fo- 
cused above on the net condensation for these systems. It 
is worth noting, however, that the oscillatory components 
of these potentials can introduce attractive terms to the 
interaction between two shifted, identical molecules, just 
as in the charge modulated example presented earlier. 



V. DISCUSSION 

In this paper, a collection of simple examples have been 
presented which demonstrate how inherent modulations 
can alter the potentials within electrolyte solutions. It 
was first shown that charge modulation can result in like- 
charge attraction between identical molecules within the 
P-B formalism. The purpose of including this discussion 
was only to make the technical point that attraction is 
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possible even when the free screening charges are consid- 
ered at the mean-field level, contrary to some suggestions 
in the literature. Next, it was demonstrated that charge 
modulation results in extra screening charge condensa- 
tion when the P-B equation is expanded to second order 
in the potential. This result, when taken together with 
the results of Ref . [T9| , implies that charge modulation- 
induced condensation universally occurs within the P-B 
formalism when the potential is expanded to second order 
in the modulation amplitudes. In particular, the result 
applies to systems which have large modulations but little 
net charge. In systems such as these, charge modulations 
can easily result in overcharging. It was further shown 
that screening charges of higher valence allow for a reduc- 
tion in the net cancellation of modulation-induced con- 
densation which occurs in neutral, univalent electrolyte 
systems. This is consistent with the fact that overcharg- 
ing has only been observed in systems which contain mul- 
tivalent counterions. Similar results were shown to apply 
for surface modulated systems and related effects can be 
expected whenever non-linearities appear in either the 
differential equation or the boundary conditions which 
determine the potential. 

We conclude by elaborating upon how one can test 
whether charge modulation significantly contributes to 
any experimental observation of overcharging. In the 
above we noted that mobile charges of equal valence but 
opposite sign are expected to condense in equal num- 
bers onto a substrate due to inherent charge modula- 
tions. This is not the case for fluctuation-induced con- 



densation, however, where charges are only expected to 
condense onto substrates of the opposite sign. This sug- 
gests that an overcharging observation can be identi- 
fied as modulation-induced if the overcharging is negated 
through the introduction of equal quantities of multiva- 
lent counterions of each sign into the system. We note 
that this test would not apply to systems where small 
surface modulations cause condensation in the weak po- 
tential limit, where the linearized P-B equation applies. 
It would, however, also apply to systems where surface 
modulations bring about condensation at higher orders 
in the P-B expansion. 



Note added to manuscript 

Following the publication of this article, I became 
aware of two references relevant to the included dis- 
cussion regarding interactions between heterogeneously 
charged surfaces [2^, [2^. Both of these predate the 
present article, and correctly indicate that charge modu- 
lation can allow for attractions between like-charged sur- 
faces. 
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